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Abstract 

We present a simple derivation of the Lamb shift using effective field theory tech- 
niques and dimensional regularisation. 
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The explanation of the Lamb shift is one of the early applications of the quantised 
electromagnetic field 0. Although many text books contain detailed expositions on this 
effect, the conventional derivations are subtle and rather intrincated. In particular a careful 
separation of small and large momenta must be carried out at some stage of the calculation. 

It has already been noticed that modern techniques of effective field theory lead to 
a simpler derivation of the Lamb shift. The important point to be realised is that in non- 
relativistic bound states there is a hierarchy of well separated scales Namely, the mass 
of the particle which forms the bound state (hard), the typical relative momentum in the 
bound state (soft) and the typical energy of the bound state (ultrasoft). In ref. counting 
rules were given to work out the size of a given NRQED diagram in the Coulomb gauge for 
soft and ultrasoft photons respectively. These rules were used in ref. to obtain the Lamb 
shift 0{ma^), which requires the matching of QED to NRQED at one loop. The calculations 
were done using a cut-off and a photon mass to regulate the UV divergencies in the effective 
theory and IR divergences in both theories respectively. 

In this note we point out that a further step in clarity and simplicity is achieved by 
(i) matching NRQED to an effective theory for ultrasoft photons which we have called 
potential NRQED (pNRQED) and (ii) using dimensional regularisation (DR) for both 
UV and IR divergences. Since we are interested in the physics at the ultrasoft scale, we may 
sequentially integrate out energies and momenta at the hard and soft scales. This is carried 
out by matching to suitable effective field theories. Integrating out the hard scale leads from 
QED to the celebrated NRQED. Integrating out the soft scale produces potential terms and 
leads from NRQED to pNRQED. In both cases the matching is most efficiently done in DR. 

The matching from QED to NRQED at one loop for the bilinear term in fermions has 
been carried out using dimensional regularisation in and for the four-fermion operators 
in 0. Since we are interested in the interaction with a static source the infinite mass limit 
for the nucleus field is taken. Hence the only relevant terms at 0{ma^) in the NRQED 
Lagrangian are the following. 
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where zD" = ido — eJ^ , zD = zV + eA on the electron field ^/^ and iD^ = ido + eZA^ on 
the nucleus field. The matching coefficients read 
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The fact that only the terms displayed in (1) are relevant can be easily seen by drawing 
diagrams of one electron one nucleus irreducible Green function (i.e. diagrams which cannot 
be disconnected by cutting one electron and one nucleus line) and taking into account that 
the next relevant scale is ma. For diagrams which cannot be disconnected by cutting a 
photon line the order of the leading contribution of each diagram is ma'""''*^^, where r is the 
number of a that appear explicitly in the diagram and s the number of 1/m factors (which 
must be compensated by ma until we obtain dimensions of energy). For diagrams which can 
be disconnected by cutting a photon line there is an extra suppression if n time derivatives 
act on this photon line. This is due to the fact that these time derivatives are only sensible 
to the typical energy. The extra suppression factor is a". We should keep in mind however 
that a given diagram in NRQED in general contains subleading contributions in a as well 
On the contrary pNRQED will be specially designed in such a way that each term and 
diagram contributes to a single order in a. 

Before going on, let us briefly comment on di in (1). This factor can be set to one by 
rescaling the photon field. The only effect this rescaling has in the effective Lagrangian is 
redefining the electric charge e. The redefined e is the low energy electron charge which 
is the one actually measured in low energy experiments, namely the one whose value is 
e^/An = a ~ 1/137. Therefore di will be put to one from now on. 

pNRQED must be regarded as an effective field theory where electron and photon energies 
and (relative) momentum ~ ma have been integrated out. It is important to realise that 
because both the energy and momentum we are integrating out are ~ ma we can use HQET 
(static) propagators for the electron field to carry out the matching. This allows to match 
NRQED and pNRQED at a given order of 1/m and a. The pNRQED Lagrangian can be 
written down at two levels: (i) in terms of electron and nucleus fields and (ii) in terms 
of bound state wave function fields. The level (i) is suitable for the matching calculation 
whereas the level (ii) is more convenient for bound state calculations. In fact it is at this level 
where the order of each term becomes explicit. The pNRQED (i) Lagrangian is a functional 
of i{j{t,x.), N{t,x.) and ultrasoft photon fields y4^(t, x), which also contains (potential) terms 
non-local in space. y4^(t,x) being ultrasoft means that they must be multipole expanded 
about the position of the nucleus in the one electron one nucleon Green function. Recall 
that this is due to the fact that in a bound state x ~ 1/ma whereas the typical energy and 
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momentum scales for ultrasoft photons are ~ ma^. 

The matching for the electron and nucleus bilinears is trivial since the related Green func- 
tions are blind to the relative momentum. Hence the terms bilinear in electron fields and the 
terms bilinear in nucleus fields in NRQED remain the same in pNRQED. The terms contain- 
ing photon fields only also remain exactly the same in NRQED and pNRQED. However we 
should keep in mind that now they correspond to ultrasoft photons only. The matching of 
one electron one nucleus Green functions induces potential terms in pNRQED. Namely terms 
bilinear in the electron and nucleus fields which are local in time but non-local in space. If we 
use the Goulomb gauge, at the order we are interested in it is enough to calculate the above 
Green functions at zero energy and at some non-zero relative momentum p. Indeed, if the 
photon fields in this Green function calculated from the pNRQED Lagrangian are multipole 
expanded, any loop gives zero in DR as we argue next. First of all there is no scale for the 
energy. Second, there is a scale for the momentum but both the HQET propagators and the 
(multipole expanded) transverse gluon propagators are insensitive to it. Hence all loops in 
pNRQED have no scale and can be put to zero in DR. Consequently the potential terms in 
pNRQED can be read off directly from the calculation in NRQED. 

A word of caution is needed. Strictly speaking the procedure above is not complete: it 
misses off-shell photons of energy ~ ma"^ and momentum ~ ma. In particular, we miss 
the contribution of diagrams which can be disconnected by cutting a photon line when time 
derivatives appear in this line. In the Coulomb gauge such photons start playing a role at 
order ma^. However, if a covariant gauge is used longitudinal photons in that region start 
playing a role already at order ma'^ and hence they must be taken into account to calculate 
the Lamb shift in those gauges. Anyway, the matching procedure can be refined so that the 
effect of such photons is included. We shall elaborate on this point elsewhere. 

At the order we are interested in only the NRQED diagrams in Fig. 1 contribute to the 
matching. We obtain 
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It is convenient to project the Lagrangian (at level (i)) above to the one electron one 
nucleus subspace and thus obtain the pNRQED Lagrangian at level (ii). This subspace is 
spanned by 

j d=^Xid=^X2^(Xi,X2)V^(Xi)Art(x2)|0) , (4) 
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where |0) is the subspace of the Fock space containing zero electrons and nuclei but an 
arbitrary number of ultrasoft photons. Since there are not spatial derivatives acting on X2 
the nucleus can be approximated by a static source which may be taken at the origin. Then 

(/?(Xi, X2) = (/?(Xi)(5^(x2) . (5) 

The dynamics of (p{t,:x.) is then described by the following Lagrangian 
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where we have multipole expanded y4^(t,xi) about zero. Each term in this Lagrangian has 
a well defined size, di acting on the electron field and l/|x| are ~ ma. Space derivatives 
acting on >l^(t,0), time derivatives and ^4^(^,0) itself are ~ ma'^. The contribution of a 
given diagram to the energy is given by ma'^'^^^'^^'^'^ . r is the number of explicit a in the 
diagram, s the number of 1/m and p the number of minus three times the number of 5^(x) 
minus the number of l/|x| minus the number of 9j. We have kept in (6) only the terms which 
are relevant to calculate the energy at 0{ma^). There are more terms of the same size (e.g. 
— ie(x.VA(i, 0)).V/m) which however do not contribute to the energy at the desired order 
as it can be seen by counting the order of the diagrams in which they appear. (Dropping 
the term —ex^xWidjAQ{t, 0) requires a more detailed analysis if a covariant gauge is to be 
used) . 

Gauge invariance can be checked order by order in this counting by introducing a gauge 
covariant field S(t,x.) as V9(t,x) = [e*^ Jo 5(i,x), such that it transforms S{t,x.) 
g-«e(z-i)0(t,o)^^^^ From this transformation it is clear that we should regard iS'(t, x) like 
an ion field rather than an electron field. After multipole expanding at the desired order, we 
obtain an explicitly gauge invariant Lagrangian 
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It is also interesting to notice that the same Lagrangian above can be obtained by using 
backwards the equation of motion in the first term of the second line in (6). 
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j d^xipH—A{t,0) ■ 'Vcp — j d'^x(pUeA{t,0) ■ x,ho (p = j d'^xip^edoA{t,0) ■ xcp , 

(8) 

where 

The Lamb shift receives contributions from the propagation of ultrasoft photons in the 
bound state. Notice that since the lagrangian of pNRQED (ii) is exphcitly gauge invariant we 
can use any gauge to calculate with it. Still the Coulomb gauge continues to be advantageous. 
In this gauge AQ{t,0) can be dropped if we use DR since it only produces tadpoles. If a 
covariant gauge is used the propagation of ^o(^, 0) must also be taken into account. In the 
Coulomb gauge at 0{ma^) the ultrasoft photons contribute to the bound state energy only 
through the pNRQED diagram of Fig. 2. In order to calculate this contribution we consider 
the two point function 

n(g,x) := I < T{ip{0)ip\x)} > (10) 

when g — > We write 
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n(g, x) must be calculated in D dimensions. By introducing the complete set of eigen- 
function oi ho in D — 4: — e dimensions ^o^n — En4>m we obtain 

nfex) ~ ':mm _ e= *m ^in\.'Hi,(, - £„)(„|v^i„)Mm , (12) 
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where v = —i^jm. We have then to calculate the following integral 
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This integral is carried out by (i) first integrating over ko and then (ii) choosing a contour 
in the |k| complex plain which embraces the cut on the positive real axes and closes itself at 
infinity in such a way that the pole of the integrand is always left inside the contour. Notice 
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that the divergent part of this integral is a polynomial. Before going on let us check that it 
can be absorbed by a counterterm in LpjvRQSD- Indeed, 

^us,iocai^^ ^ ^(n|v|m) (g - Em) (m|v|n) = (n|v (g - ho) v|n) (14) 

m 

= (^1 - ^o) + ^v^ (q -ho) [v [v, q - Kq]] ) \n) . 

In order to identify SE^ the strict limit q ^ En must be taken and hence the first two 
terms in the last expression can be dropped. The remaining term is nothing but the D — 1 
dimensional Laplacian acting on the D dimensional Coulomb potential which gives a D — 1 
dimensional (5-function. Thus (14) reads 

SUS,iocaij^^ ^ (^-|^5^-i(x)) \n) , (15) 

which can be absorbed by a counterterm in the second last term of Lp^rqed in (6) and 
(7). Since we have used MS subtraction scheme in the matching we must use the same 
subtraction scheme here. Since there are no further divergences left the limit D ^ A can 
be safely taken. Notice that this procedure avoids calculating exphcitly 0„ and En in D 
dimensions. We finally obtain that the contribution to the energy given by the diagram in 
Fig. 2 is 
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Prom the imaginary part of (13) we also obtain the total width 
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rn = E ^«l(r^|v|m)|2 {En - Em) . (17) 



The remaining contribution at 0{ma ) arises from the above mentioned potential terms 
in (6) and (7). It reads 

5' En = S^'^'En + 6^'' En + S^'^ En , 

(18) 
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(In) = \nlj)). Thus the correction to the bound state energy 6 En up to ma^ included is 
obtained by adding up the soft and ultrasoft contributions 

5£;„ = (5^E„ + 5^^E„, (19) 

which agrees with the well known result. Recall that the hard contributions are encoded 
in cd, cs and d2. Notice also that the subtraction point dependences of (in cd) and 

6^^ En cancel out, as they should. However we should stress that in order to obtain the finite 
pieces right, it is necessary that the same subtraction scheme, namely MS, which has been 
used in the matching be used in the calculations of pNRQED. 

In summary we have presented a simple and complete derivation of the Lamb shift based 
on modern EFT techniques and dimensional regular isat ion. The derivation goes through two 
EFTs, namely NRQED and pNRQED, so that at any stage it becomes clear which terms and 
diagrams are to be considered to carry out a calculation at the desired order. The formulation 
of pNRQED, namely an effective field theory for ultrasof photons, has received quite some 
attention lately [|, |. In ref. the matching from NRQED (NRQCD) to pNRQED 
(pNRQCD) was outlined and the pNRQED and pNRQCD Lagrangians for positronium 
and quarkonium were presented. Here we have worked out the pNRQED Lagrangian for 
Hydrogen-like atoms and presented its first application. Gauge invariance is manifest in 
the two EFTs. This is particularly remarkable in pNRQED, if we take into account that 
most of the previous attempts to deal with ultrasoft photons were strongly based on the 
Coulomb gauge. Last but not least DR allows to regulate both UV and IR divergences 
in a manifestly gauge invariant way (unlike the photon mass) and makes the matching 
calculations straightforward. These techniques are most promising when applied to bound 
states of heavy quarkonia in QCD [0. 
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Captionl. The first diagram is the Coulomb potential. The circle is the vertex pro- 
portional to C/), the square to Cs (spin dependent) and the dashed dot to di (the vacuum 
polarization). The dashed, solid and double lines are the static photon, electron and nucleus 
propagators respectively. 

Caption2. The thick fine and wavy fine are the ion and transverse photon propagators 
respectively. 
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Figure 1: 




Figure 2: 
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